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INTRODUCTION 


One  of  the  serious  problems  in  satellite  infrared  imagery  analyses  and  in 
the  retrieval  of  temperature  profiles  in  the  presence  of  clouds  from  sat- 
ellite soundings  is  how  to  differentiate  between  the  derived  temperature 
values  which  arise  from  the  radiation  field  involving  the  blackbody  or  non- 
blackbody  cloud.  This  problem  has  been  recognized  by  Glahn  [1]  as  one  of 
several  difficult  problems  in  conjunction  with  the  determination  of  cloud 
top  heights  and  areal  coverage.  Chahine  [2,3,4],  among  others,  in  a series 
of  outstanding  papers  has  described  methods  for  retrieving  atmospheric 
temperatures  in  the  presence  of  clouds.  However,  these  methods  assume 
that  the  clouds  are  blackbodies.  Although  a number  of  papers  [5, 6, 7, 8] 
have  examined  the  approximate  thickness  required  for  a cloud  to  radiate 
as  a blackbody,  the  cloud  blackness  problem  has  not  yet  comprehensively 
been  investigated.  If  the  cloud  thickness  is  the  required  parameter  to 
determine  whether  a cloud  temperature  arises  from  a blackbody  radiation 
field,  then  it  seems  that  the  problem  is  far  from  being  satisfactorily 
resolved. 

The  purpose  of  this  paper  is  to  discuss  the  problem  of  the  cloud  blackness 
in  the  llpm  infrared  window  region  in  relation  to  satellite  application. 

In  examining  the  internal  structure  of  the  cloud  radiation  field,  a method 
was  developed  by  which  the  degree  of  the  cloud  blackness  may  be  derived  by 
adding  a 10.5pm  channel  to,  for  example,  the  existing  12.0pm  channel  in 
the  vertical  temperature  profile  radiometer  (VTPR).  Consequently,  it 
would  be  possible  to  determine  the  reliability  of  the  derived  temperature 
values.  To  develop  the  method,  it  is  necessary  to  first  understand  phys- 
ically the  internal  cloud  radiation  field.  Since  the  primary  interest  was 
in  satellite  application,  only  the  vertical  upward-going  radiation  of 
clouds  was  investigated. 

The  computer  program  for  solving  the  radiative  transfer  equation  was  based 
upon  the  Gauss-Seidel  iteration  technique  [9,10]  and  had  been  checked  out 
against  published  values  in  the  literature  [5,8].  For  water  vapor  attenu- 
ation, use  was  made  of  LOWTRAN  3 [11]  and  the  empirical  formula  given  by 
Roberts  et  al.,  [12]  to  account  for  selective  and  continuum  absorption, 
respectively.  The  computer  program  not  only  gives  the  cloud  radiance  values 
as  a function  of  zenith  angles  at  different  prescribed  geometric  levels, 
but  also  calculates  the  optical  properties,  including  the  optical  thick- 
nesses and  single-scattering  albedos  at  corresponding  levels. 

Since  an  understanding  of  the  internal  structure  of  the  cloud  radiation 
field  is  essential  to  the  development  of  the  technique,  the  next  section 
describes  the  cloud  models  to  be  used  in  the  investigation.  The  theoret- 
ical basis  then  follows.  In  section  4,  computational  results  are  presented, 
and  the  internal  radiation  structures  are  examined.  These  analyses  lead 
to  the  development  of  the  technique  for  determining  the  degree  of  cloud 
"blackness"  or  emissivity  in  section  5.  A few  numerical  examples  will  be 
given  here.  The  final  section  contains  the  findings  and  conclusions. 
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CLOUD  MODEL 


Only  recently  have  the  radiative  properties  of  clouds  been  studied  realis- 
tically with  inhomogeneous  droplet  and  temperature  distributions  considered. 
Stephens  [13],  using  the  meteorological,  microphysical,  and  radiation  data 
furnished  to  him  by  Paltridge  [14],  investigated  the  inhomogeneous  strati- 
form clouds  theoretically  and  found  the  calculated  downward  radiances  in 
good  agreement  with  the  observed  values.  Platt's  observations  [15]  of  up- 
welling  radiances  emerging  at  different  levels  from  his  inhomogeneous 
Stratocumulus  Deck  2 as  shown  in  Fig.  1 were  adopted  in  this  paper.  This 
cloud  had  a thickness  of  abou^  550  m with  a temperature  distribution  as 
given  in  his  figure.  From  the  typical  humidity  profile  shown,  it  may  be 
inferred  that  except  for  the  middle  region  of  the  cloud,  which  was  satu- 
rated, both  the  top  and  base  layers  were  most  likely  subsaturated. 

Platt's  radiance  profile  in  the  10.5um  window  was  simulated  and  radiance 
values  were  expressed  in  mW  m-2  sr“i  cm,  which  can  be  readily  converted 
to  his  physical  unit  by  means  of  a simple  formula.  However,  it  should  be 
noted  that  the  choice  of  cloud  microphysics  and  wavelength  would  in  no  way 
affect  either  the  analysis  of  the  internal  radiation  structure  or  the  find- 
ings and  conclusions. 

Although  no  microphysics  data  were  provided  by  Platt  [15],  on  the  basis 
of  a number  of  cloud  physics  studies  [16,17,18],  it  appears  that  the 
gamma  distribution  [19]  with  two  parameters,  a = 4 and  8=1,  may  ade- 
quately represent  the  stratocumulus  cloud  over  the  ocean.  The  droplet 
sizes  ranged  from  0.5ym  to  25um  with  a mean  radius  of  5tjm.  In  view  of  the 
abundant  supply  of  sea-salt  nuclei  over  the  ocean,  the  upper  limit  of  the 
range  seems  to  represent  a reasonable  value.  With  the  known  distribution 
function  and  size  range  and  the  refractive  indices  for  water  spheres  [21], 
Mie  scattering  and  absorption  coefficients  per  unit  length  may  be  calcu- 
lated. 

Stratocumulus  Deck  2 was  divided  into  22  layers  each  25  m thick,  and  radi- 
ance values  were  computed  at  each  level.  With  the  Mie  coefficients  known, 
a first  approximation  to  the  number  density  of  droplets  required  to  emit 
the  desired  amount  of  radiation  at  each  level  was  made.  Then  the  tempera- 
ture profile  and  the  estimated  droplet  number  densities  at  chosen  levels 
were  fed  into  the  computer  program  to  generate  a first  estimate  of  upward 
radiances  as  a function  of  cloud  depths.  Adjustments  were  made  of  the 
number  densities  at  several  levels  in  each  run.  Platt's  radiance  profile 
was  reproduced  almost  exactly  for  che  10.5Mm  window  in  four  trial  runs. 

The  input  meteorological  and  microphysical  values  are  given  in  Table  1. 
Although  the  computer  program  is  capable  of  accounting  for  water  vapor 
attenuation  both  inside  and  outside  the  cloud,  the  attenuation  of  surface 
radiation  below  the  cloud  base  was  not  considered  since  the  amount  of 
attenuation  was  quite  small,  as  is  evident  in  Fig.  1. 


AIR  TEMPERATURE 


Figure  1.  Profiles  of  vertical  upward  radiance,  air  temperature,  and  blackbody 
radiance  through  the  stratocumulus  cloud  deck. 
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TABLE  1.  INPUT  METEOROLOGICAL  AND  MICROPHYSICAL  DATA  OF  THE  MODEL  CLOUD  I 


Height 

(km) 

Pressure 

(mb) 

Ai rtemp 
(°K) 

Dewtemp 

(°K) 

Droplet 
( cm" 3 ) 

2.25 

777.21 

280.55 

279.15 

5.0 

2.10 

791.40 

281.05 

280.65 

10.0 

2.00 

801.00 

281,65 

281.65 

15.0 

1.90 

810.48 

282.24 

282.25 

25.0 

1.85 

815.26 

282.45 

282.45 

50.0 

1.80 

820.07 

282.45 

282.35 

100.0 

1.75 

824.91 

282.35 

282.25 

130.0 

1.70 

829.78 

281.95 

281.75 

120.0 

Since  most  of  the  droplets  concentrate  in  the  lower  150  m layer,  rapid 
attenuation  of  the  surface  radiation  is  found  in  approximately  the  first 
200  m above  the  cloud  base  (Fig.  1).  The  total  number  of  droplets  in  the 
550  m column  per  unit  cross-sectional  area  is  2.20  x 10^.  These  droplets 
were  redistributed  uniformly  throughout  the  entire  column,  resulting  in 
40  droplets  cm"^  or  3.52  x 10"^  g m-^  in  liquid  water  content  at  each 
level.  By  taking  the  average  temperature  in  the  vertical  column,  a 
homogeneous  and  isothermal  cloud  similar  to  those  used  by  Yamamoto  et  al. 
[5]  and  Hunt  [8]  was  created  from  the  original  cloud  model. 


THEORETICAL  BACKGROUND 


The  equation  governing  the  transfer  of  thermal  radiation  through  a scat- 
tering and  absorptive  medium  is 


dI(T;tj)  _ 


w (t)  -+1 

— 2^ — j p°(u,m' )I(t,u' )dp' 


- [1  - w^(t)]B[T^(t)]  , (1) 

where  I is  the  spectral  intensity  or  radiance,  y the  cosine  of  the  zenith 

angle,  t the  optical  thickness,  the  single-scattering  albedo,  p (y,y') 

the  integrated  phase  function  over  azimuth  angle,  and  B[T^(t)]  the  Planck 

function  for  cloud  temperature  T^  at  level  t.  For  simplicity,  the 

wavenumber-dependent  subscript  has  been  eliminated.  The  above  equation  can 
be  separated  into  two  components  as  follows:  one  component  denotes  the 
external  field,  primarily  due  to  surface  radiation,  and  the  other  denotes 
the  cloud  emission. 
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For  the  external  field  (Ig)» 


M 


(T.y) 


= Ig(T;u) 


I 


p°(ii,M'  )lg(T,W'  )du' 


» 


with  boundary  conditions 


lg{0;-vi)  = 0 
Ig(xb;+u)  = Bd^)  ) . 


where  0 and  denote  the  cloud  top  and  base,  respectively,  and  T, 
surface  temperature. 

For  the  cloud  emission  (1^,). 


dl  (t;h) 


(y .y'  )Ij.(T;y ' )dy ' 


- [1  - 2:^(t)]B[T^(t)]  , 

with  boundary  conditions 


l^(0;-y)  =0  j ; 

^c(^b’''*^^  = 0 j . 


To  simplify  the  expressions,  let 


S;(T;y)  = -j- 


r 


p°(y,y' )Ig(T;y' )dy' 


(2) 


(3) 


is  the 


(4) 


(5) 


(6) 
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and  : 

w (t)  -+1  „ 

S'(t;u)  = -25—  I P (i'.i'')Ic(T!p')'iu'  . (7) 

^*-1 

Equations  (6)  and  (7)  represent  the  contributions  by  droplet  scatterings.  I 

Thus,  equations  for  upward  radiation  are  | 

dl  (f.+y)  I 

y — ^ j Y — = Ig(T;+y)  - Sg(T;+y)  , (8)  | 

and  ; 


dl  (x-.+y) 

u — = I(.(T;+y)  ■■  S^(T;+y) 


[1  - 2;,(x)]B[T^(t)]. 


(9) 


[I 


■i 

I 
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corresponding  to  (2)  and  (4),  respectively.  Similar  expressions  can  be 
written  for  downward  radiation  by  changing  +u  to  -y. 

For  a finite  optical  thickness,  the  expressions  are  readily  solved  [20]. 

The  computer  program  solves  (8)  and  (9)  in  essentially  the  same  manner  by 
means  of  the  Guass-Seidel  iteration  technique  for  upweiling  (as  well  as 
down-welling)  radiation  at  successive  levels  from  the  cloud  base  to  the 
top  in  accordance  with  boundary  conditions  (3)  and  (5). 

The  approximate  solutions  of  (8)  and  (9)  with  the  boundary  conditions 
specified  for  the  upweiling  radiances,  Ig  and  1^.,  at  each  successive 

level  from  the  cloud  base  may  be  written  in  the  form 

1 


I 
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r B T + S 
9 


9 


(10)  j 

i 


Ml  - - T)  + , 


(11) 


where  represents  B(Tg);  B^  B[T^(t)];  Sg  S^(l  - T);  S^(l  - T); 

and  T is  the  vertical  transmittance  from  the  cloud  base  to  that  level. 

Values  of  Ig  and  for  the  cloud  model  were  obtained  for  lO.Sym,  ll.Oym, 

11.5ym,  12.0ym,  and  12.5ym.  The  refractive  indices  for  these  wavelengths 
were  taken  from  Irvine  and  Pollack  [21].  The  volume  scattering,  absorption. 
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and  extinction  coefficients  (per  km)  for  100  droplets  and  the  average 
single-scattering  albedo  for  the  model  cloud  are  given  in  Table  2. 


Upwelling  radiances  1^  and  together  contribute  to  the  total  upwelling 
radiation  from  the  model  cloud, ’whereas  represents  radiance  from  an 
isolated  cloud  only,  i.e.,  no  surface  emission  contribution  [61. 

TABLE  2.  VOLUME  SCATTERING  (35^,3^).  ABSORPTION  AND 

EXTINCTION  COEFFICIENTS  AND  AVERAGE  SINGLE 

SCATTERING  ALBEDO  OF  THE  MODEL  CLOUD  AS  A 
FUNCTION  OF  WAVELENGTH  (X) 


X 

(ym) 

^scat 

(km"^) 

^abs 
(km-^ ) 

®ext 

(km"^) 

'h 

^0 

10.5 

7.2609 

6.1013 

13.3622 

0.5044 

11.0 

4.8780 

7.2201 

12.0981 

0.3665 

11.5 

4.6519 

8.5988 

13.2507 

0.3174 

12.0 

5.2469 

9.7333 

14.9802 

0.3152 

12.5 

5.8936 

10.3532 

16.2468 

0.2796 

INTERNAL  RADIATION  STRUCTURE  AND  DISCUSSIONS 

Examination  of  (10)  and  (11)  reveals  that  as  T becomes  closer  to  unity 
(or  T approaches  0)  near  the  base  of  a thick  cloud  the  total  upward 
radiance,  I + I , originates  mainly  from  surface  radiation  and  little 
9 ^ 

from  cloud  emission.  On  the  other  hand,  near  the  cloud  top  where  T ap- 
proaches 0 (or  T approaches  infinity),  the  cloud  emission  and  scattering 
dominate. 

Thus, 

Ig  - 0,  then  I^  = B[T^(0)]  = B^ 
and 


6c  = <’  - “=)6c  * 6,  . (12) 
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This  means  that  the  apparent  blackbody  radiation  at  the  top  of  a "black" 

cloud  is  composed  of  two  radiant  components,  which  include  partly  thermal 
emission  and  partly  scattering.  The  fractional  amount  of  each  is  deter- 
mined by  the  average  single  scattering  albedo  of  the  cloud. 

Since  the  values  of  all  the  terms  at  each  level  in  (10)  and  (11),  except 

those  of  S„  and  which  also  include  some  contributions  from  downward 
9 c 

radiation,  can  be  readily  calculated,  it  is  simple  to  sort  out  the  contri- 
butions to  the  vertical  upward  radiation  made  by  cloud  emission,  surface 
radiation,  and  droplet  scatterings  and  to  examine  their  changes  with  wave- 
length and  cloud  depth.  For  simplicity  of  presentation,  the  internal 
cloud  radiation  field  was  examined  at  100  m intervals  up  to  500  m as  well 
as  the  cloud  top  at  550  m.  Also,  the  cloud  model  was  examined  in  four 
different  cases  involving  (1)  inhomogeneous  droplet  and  temperature  dis- 
tributions, (2)  homogeneous  in  both  droplet  and  temperature  distributions, 
(3)  isothermal  but  inhomogeneous  in  droplet  distribution,  and  (4)  homoge- 
neous in  droplet  distribution  but  with  nonisothermal  temperature  distri- 
bution. The  first  case  was  run  for  all  the  wavelengths  listed  in  the  pre- 
ceding section,  and  the  other  three  at  10,5ym  and  12,5ym  only.  In  this 
manner,  some  physical  insight  may  be  derived  as  to  the  suitability  of 
employing  a homogeneous  cloud  as  opposed  to  an  inhomogeneous  one  for  mod- 
eling purposes. 

Three  tables  are  presented.  The  first  two  deal  with  percentage  contribu- 
tions in  the  four  cases,  in  which  SS  denotes  the  total  scattering  contri- 
butions, i,e,,, Sg  + S^,  CE  the  cloud  emission,  and  TS  the  reduced  surface 

radiation.  To  evaluate  what  role  the  inclusion  of  water  vapor  would  play 
in  cloud  modeling,  the  calculated  radiance  values  at  each  level  are  given 
in  Table  5,  so  that  the  percentage  contribution  made  by  water  vapor  to 
upward  radiation  may  be  understood. 

In  Table  3,  the  changes  with  wavelength  of  percentage  values,  when  compared 
with  those  of  optical  values  in  Table  2,  show  that  (1)  the  surface  contri- 
butions (TS)  at  each  level  vary  inversely  with  the  extinction  coefficients 
(Bext^’  (2)  the  total  scattering  contributions  (SS)  fluctuate  about  the 

average  single  scattering  albedos  (w^)  of  the  cloud,  and  (3)  the  emission 

contributions  (CE)  increase  with  the  wavelength.  An  examination  of  Tables 
3 and  4 shows  that,  at  all  wavelengths,  as  the  surface  contributions 
decrease  with  thickness,  the  emission  contributions  increase.  The  scat- 
tering contributions  due  to  surface  radiation  and  those  due  to  cloud  emis- 
sion also  behave  in  similar  manners,  although  they  are  not  tabulated  here. 
However,  the  percentage  values  of  the  total  scattering  contributions  in 
the  four  cases  at  corresponding  wavelengths  are  vastly  different,  except 
for  case  3,  This  can  be  explained  by  the  fact  that  cloud  microphysical 
properties  play  significant  roles  in  the  determination  of  the  internal 
radiation  structures. 
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TABLE  3.  PERCENTAGE  CONTRIBUTIONS  BY  TOTAL  SCATTERING  (SS), 
CLOUD  EMISSION  (CE),  AND  REDUCED  SURFACE  RADIATION 
(TS)  TO  UPWELLIN6  RADIATION  AT  DIFFERENT  LEVELS  (tn) 
AS  A FUNCTION  OF  WAVELENGTH  IN  CASE  1. 


Wavelength(iJtn) 

Level 

(m) 

10.5 

11.0 

11.5 

12.0 

12.5 

100 

SS 

44 

32 

29 

32 

31 

CE 

34 

43 

49 

52 

54 

TS 

22 

25 

22 

18 

15 

200 

SS 

50 

36 

32 

32 

32 

CE 

40 

51 

57 

60 

62 

TS 

10 

13 

11 

8 

6 

300 

SS 

50 

36 

32 

32 

31 

CE 

42 

54 

60 

62 

65 

TS 

8 

10 

8 

6 

4 

400 

SS 

50 

36 

32 

32 

30 

CE 

43 

55 

61 

63 

67 

TS 

7 

9 

7 

5 

3 

500 

SS 

50 

36 

31 

31 

28 

CE 

44 

56 

63 

65 

70 

TS 

6 

8 

6 

4 

2 

550 

SS 

49  ‘ 

36 

31 

31 

28 

CE 

45 

58 

64 

65 

70 

TS 

6 

6 

5 

4 

2 
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TABLE  4.  PERCENTAGE  CONTRIBUTIONS  BY  TOTAL  SCATTERING,  CLOUD 
EMISSION,  AND  REDUCED  SURFACE  RADIATION  TO  UPWELLING 
RADIATION  AT  DIFFERENT  LEVELS  AT  lO.Sum  AND  IZ.Sum 
IN  CASE  2 (HOMO  DT),  CASE  3 (HOMO  T ONLY),  AND  CASE 
4 (HOMO  D ONLY). 


Level  (tn) 

Case  2 

Homo  DT(pm) 
10.50  12.50 

Case  3 

Homo  T Only(pm) 
10.50  12.50 

Case  4 

Homo  D Only(iim) 

10.50  12.50 

100 

SS 

23 

17 

44 

31 

23 

17 

CE 

17 

32 

34 

54 

17 

32 

TS 

60 

51 

22 

15 

60 

51 

200 

SS 

36 

25 

50 

32 

36 

25 

CE 

28 

49 

39 

62 

28 

49 

TS 

36 

26 

11 

6 

36 

26 

300 

SS 

34 

29 

51 

31 

34 

29 

CE 

29 

58 

41 

65 

30 

59 

TS 

37 

13 

8 

4 

36 

12 

400 

SS 

48 

31 

50 

30 

48 

31 

CE 

39 

63 

43 

67 

39 

63 

TS 

13 

6 

7 

3 

13 

6 

500 

SS 

51 

31 

50 

28 

50 

31 

CE 

42 

66 

44 

70 

42 

66 

TS 

7 

3 

6 

2 

8 

3 

550 

SS 

51 

32 

49 

28 

51 

31 

CE 

43 

66 

45 

70 

43 

67 

TS 

6 

2 

6 

2 

6 

2 
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TABLE  5.  PERCENTAGE  EFFECT  OF  WATER  VAPOR  ON  UPWELLING  RADIATION 
AT  DIFFERENT  LEVELS  AS  A FUNCTION  OF  WAVELENGTH. 


Level 

10.50 

11.00 

11.50 

12.00 

12.50 

(m) 

Vapor 

A(iJin) 

A(Min) 

A()jm) 

A(ym) 

X(ym) 

100 

Yes 

89.34 

100.9 

105.8 

110.2 

No 

89.49 

95.77 

101.1 

106.0 

110.7 

Percentage 

-0.2 

-0.2 

-0.2 

-0.2 

-0.4 

200 

Yes 

86.48 

92.76 

98.14 

107.6 

No 

86.67 

92.96 

98.33 

103.3 

108.0 

Percentage 

-0.2 

-0.2 

-0.2 

-0.2 

-0.4 

300 

Yes 

85.52 

91.84 

97.27 

102.3 

106.8 

No 

85.78 

92.11 

97.51 

102.5 

107.2 

Percentage 

-0.3 

-0.3 

-0.2 

-0.2 

-0.4 

400 

Yes 

84.88 

91.21 

96.66 

101.7 

106.1 

No 

85.20 

91.55 

96.95 

102.0 

106.7 

Percentage 

-0.4 

-0.4 

-0.3 

-0.3 

-0.6 

500 

Yes 

84.41 

90.75 

96.20 

101.2 

105.5 

No 

84.79 

91.15 

96.56 

101.6 

106.3 

Percentage 

-0.4 

-0.4 

-0.4 

-0.4 

-0.8 

550 

Yes 

84.23 

90.58 

96.03 

101.1 

105.3 

No 

84.64 

91.01 

96.42 

101.4 

106.1 

Percentage 

-0.5 

-0.5 

-0.4 

-0.4 

-0.8 

.4 


13 


Percentage  contributions  are  almost  identical  at  corresponding  levels  at 
both  the  10.5vim  and  IZ.Spm  wavelengths  in  cases  1 and  3 which  share  the 
same  inhomogeneous  droplet  distribution.  The  percentage  values  are  nearly 
identical  in  cases  2 and  4 which  share  the  same  homogeneous  droplet  distri- 
bution. Apparently,  this  indicates  that  the  upward  radiation  emerging  from 
any  level  inside  as  well  as  at  the  top  of  a cloud  is  virtually  independent 
of  its  internal  temperature  structure;  moreover,  this  implies  that  the  so- 
called  cloud-top  temperature  may  be  regarded  as  the  mean  temperature  of  the 
cloud. 

When  the  cloud  is  viewed  as  a whole,  i.e.,  550  m in  thickness,  the  amounts 
in  all  four  cases  of  upward  radiation  emerging  from  the  top  are  nearly  the 
same  at  their  corresponding  wavelengths.  The  calculated  radiance  values 
are  84.23,  83.85,  83.95,  and  83.67  mW  m"^  sr~^  cm  at  10.5ym  for  the  four 
cases,  respectively.  At  12.5pm,  they  are  105.30,  105.50,  105.70,  and 
104.85,  respectively.  Thus  considering  the  cloud  as  a whole,  little  dif- 
ference exists  at  a given  wavelength  in  upward  radiation  emerging  from 
its  top  between  a homogeneous  model  and  an  inhomogeneous  one  so  long  as 
they  have  the  same  total  droplet  number  or  the  same  total  liquid  water 
amount  and  the  same  mean  temperature  in  the  cloud  column. 

In  the  present  study,  since  surface  emission  (TS)  still  contributes  to  the 
upward  radiance  at  the  cloud  top,  the  cloud  is  not  "black."  It  is  even 
less  black  at  10.5pm  than  at  12.5pm.  The  calculated  average  cloud  verti- 
cal emissivity  over  the  10.5pm  to  12.5pm  interval  is  0.85  compared  with 
Platt's  value  [15]  of  about  0.90  in  the  10pm  to  12pm  band.  Furthermore, 
it  may  be  observed  in  Table  3 that  at  the  550  m level  the  scattering 
contributions  (SS)  at  all  wavelengths  nearly  equal  the  average  single 
scattering  albedos  of  the  cloud  in  Table  2 at  the  corresponding  wave- 
lengths. If  the  cloud  were  to  grow  another  200  to  300  m thicker,  it 
could  be  expected,  as  indicated  in  (12),  that  at  10.5pm,  for  example, 
the  upward  emission  and  scattering  contributions  would  split  about  50-50. 
And  this  is  when  cloud  emissivity  would  reach  about  unity. 

Table  5 shows  the  effect  of  water  vapor  on  upward  radiation  at  different 
levels  as  a function  of  wavelength.  Although  there  was  appreciable  water 
vapor  inside  the  cloud  with  about  98%  relative  humidity  on  the  average 
as  given  in  Table  1,  the  presence  of  water  vapor  appears  to  have  negative 
but  negligible  effects  on  the  upward  radiance. 


RATIO  TEST  OF  CLOUD  BLACKNESS 

Not  until  a cloud  behaves  like  a blackbody  radiator  can  the  study  of  its 
internal  radiation  field  ignore  the  radiation  contribution  from  the  under- 
lying surface.  The  surface  temperature  is,  in  general,  appreciably  higher 
than  the  cloud  temperature,  while  the  cloud-top  temperature  is  lower  than 
the  base  temperature.  Only  when  the  cloud  becomes  "black"  will  the  up- 
welling  radiation  directly  relate  to  the  cloud-top  temperature.  For  a 
nonblack  cloud,  the  cloud  radiation  field  is  contaminated  by  the  surface 
radiation;  one  might  say  that  the  radiation  field  is  distorted.  To  illus- 
trate such  distortion,  Table  6 shows  the  equivalent  blackbody  temperatures 
of  the  radiation  field  at  different  levels  as  a function  of  wavelength. 
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TABLE  6.  EQUIVALENT  BLACKBODY  TEMPERATURES  (°C)  CORRESPONDING  TO 

UPWELLING  RADIANCES  AS  A FUNCTION  OF  WAVELENGTH  AT  DIFFERENT 
LEVELS  FROM  THE  CLOUD  BASE  WITH  MEAN  TEMPERATURES  (T  °C). 


Level 

(m) 

T 

(°C) 

10.50 

(pm) 

11.00 

(pm) 

11.50 

(pm) 

■t— 

1 

100 

9.13 

15.03 

14.34 

13.50 

12.88 

12.37 

200 

9.19 

13.08 

12.47 

11.72 

11.15 

10.72 

300 

9.06 

12.42 

11.85 

11.15 

10.63 

10.20 

400 

8.84 

11.98 

11.43 

10.75 

10.24 

9.75 

500 

8.62 

11.65 

11.12 

10.44 

9.92 

9.37 

550 

8.52 

11.53 

11.00 

10.33 

9.85 

9.24 

Even  though  the  mean  temperature  at  a given  level  is  at  a constant  temper- 
ature, the  equivalent  blackbody  temperatures  decrease  with  wavelength. 

Normally,  blackbody  radiance  values  increase  monotonically  with  wavelength 
at  a given  temperature.  In  the  present  situation,  the  radiance  values  at 
a given  level  in  Table  5 also  increase;  however,  the  amount  of  increase  is 
relatively  small.  The  percentage  increases  in  blackbody  radiance  at,  for 
instance,  9‘’C  are  about  9%,  17%,  24%,  and  30%,  and  at  14°C  about  8%,  16%, 

22%,  and  28%,  respectively,  for  ll.Oym,  ll.Sum,  12.0um,  and  12.5pm,  with 
reference  to  the  blackbody  value  for  10.5pm.  By  comparison,  the  percentage 
increases  in  Table  5 at,  for  instance,  the  lOO-m  level  are  about  7%,  13%, 

18%,  and  23%,  respectively.  At  the  550-m  level,  they  are  about  8%,  14%, 

20%,  and  25%,  respectively.  It  is  this  distortion  of  the  temperature-  ^ 

wavelength  relationship,  or  rather  of  the  radiance-wavelength  relationship  j 

inside  the  cloud,  that  forms  the  basis  upon  which  the  proposed  method  for  j 

determining  the  quality  of  cloud  blackness  or,  to  some  extent,  cloud  1 

emissivity  is  developed.  | 

j 

Tables  3 and  6 show  an  apparent  excess  in  the  amount  of  vertical  upward  | 

radiation,  more  at  lO.Spm  than  at  any  other  wavelengths,  giving  rise  to  I 

higher  equivalent  blackbody  temperatures  in  Table  6 at  all  levels.  On  the  \ 

other  hand.  Table  3 appears  to  indicate  that  the  surface  (TS)  and  hence  j 

the  resulting  scattering  contributions  (SS)  together  are  somewhat  greater 
at  10.5pm. 

i 

The  former  is  governed  by  the  extinction  coefficient  of  the  cloud  and  the  ' 

latter  by  its  single  scattering  albedo.  If  the  ratio  of  the  observed 

radiance  at  10.5pm  to  that  at  a longer  wavelength  in  the  same  window 

region  is  taken,  this  ratio  will  enhance  the  presence  of  an  excessive 

amount  at  10.5pm  of  surface  and  scattered  radiation  when  the  cloud  is 

not  "black."  However,  when  the  cloud  is  "black,"  i.e.,  no  contribution  ! 

from  surface  and  scattered  radiation,  the  ratio  would  be  equivalent  to 
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that  of  blackbody  radiances.  Table  7 gives  the  ratios  of  blackbody  radi- 
ances as  a function  of  temperature.  The  ll.Oum  wavelength  was  not  used, 
not  only  because  it  lies  too  close  to  the  10,5ym  wavelength,  but  also 
because  the  advantage  of  its  lower  single  scattering  albedo  is  somewhat 
offset  by  the  disadvantage  of  its  lower  extinction  coefficient.  On  the 
basis  of  these  discussions,  the  two  criteria  for  choosing  the  shorter 
wavelength  are:  (1)  relatively  higher  single  scattering  albedo  and  (2) 
relatively  lower  extinction  coefficient.  The  ll.Oym  wavelength  satisfies 
(2)  but  not  quite  (1),  as  indicated  in  Table  2, 


TABLE  7.  RATIOS  OF  BLACKBODY  RADIANCES  AT  TWO  WAVELENGTHS  AS  A 
FUNCTION  OF  AIR  TEMPERATURE  (*0) 


Temperature 

(°C) 

Wavelength  Ratio 

10.5/11.5 

10.5/12^.0 

10.5/12.5 

-25.0 

0.8108 

0.7455 

0.6934 

0.8184 

0.7555 

0.7053 

-15.0 

0.8257 

0.7652 

0.7168 

-10.0 

0.8327 

0.7746 

0.7281 

-5.0 

0.8396 

0.7838 

0.7391 

0.0 

0.8462 

0.7927 

0.7499 

5.0 

0.8527 

0.8014 

0.7604 

10.0 

0.8589 

0.8098 

0.7706 

15.0 

0.8650 

0.8180 

0.7806 

20.0 

0.8709 

0.8260 

0.7903 

25.0 

0.8766 

0.8337 

0.7998 

30.0 

0.8821 

0.8413 

0.8090 

35.0 

0.8875 

0.8487 

0.8181 

When  the  cloud  in  question  is  "black"  in  the  window  region,  the  ratio  of 
the  observed  radiances  will  give  the  cloud  temperature.  Moreover,  only 
under  this  condition  will  the  radiance  values  at  the  two  wavelengths,  upon 
conversion,  also  yield  nearly  the  same  temperature.  The  ratio  of  black- 
body radiances  in  Table  7 changes  rather  slowly  with  temperature  within 
narrow  limits,  and  it  does  not  take  many  extra  surface  and  scattering  con- 
tributions to  cause  the  ratio  to  deviate  appreciably  from  its  critical 
blackbody  value.  The  critical  ratio  is  unknown;  but  since  the  primary 
interest  lies  in  establishing  a reliability  criterion  for  accepting  or 
rejecting  the  temperature  derived  from  satellite  data,  the  problem  is 
to  find  a means  to  distinguish  between  black  and  nonblack  clouds.  For 
this  purpose,  a reference  ratio  instead  of  the  critical  ratio  is  sought. 

Besides  the  present  cloud  model,  a number  of  homogeneous  models  such  as 
those  used  by  Yamamoto  et  al.  [5]  and  Hunt  [8]  were  also  investigated. 
They  were  1 km  thick.  After  repeated  numerical  experiments,  it  was  found 
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that  the  ratio  of  the  blackbody  radiances  derived  from  the  surface  temper- 
ature served  quite  well  as  a reference  ratiOc  When  the  ratio  of  the 
observed  (or  more  appropriately,  synthetic)  radiances  at  a certain  cloud 
level  became  smaller  than  the  reference  ratio,  the  calculated  vertical 
emissivity  at  that  level  in  every  case  showed  better  than  0.95  in  the 
window  region.  The  surface  temperature  need  not  be  accurate.  Since  most 
clouds,  except  those  low-hanging  stratus  and  fogs,  are  generally  1 km  or 
more  above  the  ground,  a temperature  10  degrees  less  than  the  estimated 
surface  temperature  could  be  use^^  Over  the  continent,  even  lower  temper- 
atures could  be  used.  This  approach  was  tested  for  10  degrees  less;  and 
an  improvement  of  about  l%-2%  was  obtained. 

To  illustrate  the  technique  just  described  for  determining  the  degree  of 
cloud  blackness  or  cloud  emissivity.  Tables  8 and  9 may  serve  as  examples. 
Table  8 shows  the  inhomogeneous  model,  and  Table  9 one  of  Yamamoto's  [5] 
homogeneous  clouds.  This  particular  cloud  given  here  had  a surface  temper- 
ature of  30°C  and  a cloud  temperature  of  -30°C. 

In  the  case  of  Platt's  [15]  cloud  in  Table  8,  the  cloud  never  becomes 
"black"  since  the  observed  ratios  at  any  levels  never  become  smaller  than 
their  respective  reference  ratios,  although  at  550  m the  former  nearly 
approach  the  latter.  In  Table  9 the  ratios  are  tabulated  up  to  550  m as 
well  as  at  1 km,  the  top  of  the  model  cloud.  The  reference  ratios  at 
30°C  can  also  be  directly  obtained  from  Table  7.  The  "observed"  ratios 
at  the  cloud  top,  when  compared  with  the  corresponding  values  in  Table  7, 

indicate  that  the  equivalent  cloud  temperature  is  below  -25°C,  which  is 
compatible  with  the  given  cloud  temperature  of  the  model.  These  ratios 
are  appreciably  smaller  than  the  reference  ratios.  The  small  arrows  in 
Table  9 show  where  the  "observed"  ratios  become  smaller.  It  is  evident 
that  they  occur  at  different  levels,  450  m,  500  m,  and  550  m.  Past  studies 
[6]  show  that  cloud  emissivity  (or  blackness)  is  frequency-dependent. 

Hence,  a nonblack  cloud  at  10.5iim  may  appear  "black"  when  viewed  at  12.0um 
or  12.5um,  which  means  that  if  cloud  emissivity  is  required  to  be  the  same 
in  this  window  region  then  the  required  cloud  thickness  will  decrease  as 
wavelength  increases.  To  further  illuminate  this  point,  the  calculated 
vertical  emissivities  at  these  three  levels  are  listed  as  a function  of 
wavelength: 


Level  Wavelenqth(ym) 


(m) 

10.5 

11.0 

11.5 

12.0 

12.5 

450 

0.90 

0.93 

0.96 

0.97 

0.97 

0.92 

0.95 

0.97 

0.98 

0.98 

550 

0.94 

0.96 

0.98 

0.99 

0.99 

F 
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TABLE  8.  RATIOS  OF  THE  "OBSERVED"  RADIANCES  AT  DIFFERENT  LEVELS 
OF  THE  INHOMOGENEOUS  CLOUD  MODEL  (SURFACE  TEMPERATURE  = 
21.70*0 


Level 

(m) 

Wavelength  Ratio 

10.5/11.5 

10.5/12.0 

10.5/12.5 

550 

0.8771 

0.8331 

0.7999 

500 

0.8774 

0.8341 

0.8001 

450 

0.8778 

0.8338 

0.7999 

400 

0.8781 

0.8346 

0.8000 

350 

0.8787 

0.8351 

0.8006 

300 

0.8792 

0.8360 

0.8007 

250 

0.8801 

0.8369 

0.8018 

200 

0.8812 

0.8388 

0.8037 

150 

0.8828 

0.8408 

0.8408 

100 

0.8854 

0.8444 

0.8107 

50 

0.8848 

0.8440 

0.8110 

Surface 

0.8728 

0.8285 

0.7935 

Note : The 

last  row  contains  the 

reference  ratios. 

TABLE  9. 

RATIOS  OF  THE  "OBSERVED"  RADIANCES  AT  DIFFERENT  LEVELS 
OF  THE  HOMOGENEOUS  CLOUD  MODEL  (SURFACE  TEMPERATURE  = 
30°C  AND  CLOUD  TEMPERATURE  = -30*C) 

Level 

Wavelength  Ratio 

(m) 

10.5/11.5 

10.5/12.0 

Id. 5/12. 5 

1000 

0.8101 

0.7425 

0.6878 

550 

^0.8726 

0.8134 

0.7585 

500 

0.8895 

^0.8344 

0.7802 

450 

0.9088 

0.8593 

->0.8066 

400 

0.9297 

0.8877 

0.8376 

350 

0.9511 

0.9181 

0.8719 

300 

0.9706 

0.9481 

0.9072 

250 

0.9835 

0.9729 

0.9385 

200 

0.9913 

0.9865 

0.9586 

150 

0.9852 

0.9825 

0.9596 

100 

0.9648 

0.9565 

0.9345 

50 

0.9299 

0.9083 

0.8820 

Surface 

0.8821 

0.8413 

0.8090 

Note:  Small  arrow  points  to  where  the  observed  ratio 
becomes  less  than  the  reference  ratio. 


The  last  row  contains  the  reference  ratios. 
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If  straight  averages  are  taken  on  the  emissivity  values,  for  the  450  m, 
500  m,  and  550  m levels,  the  emissivity  values  are  0.95,  0.96,  and  0.97, 
respectively.  At  550  m Yamamoto's  cloud  appeared  "blacker"  than  Platt's. 
Although  the  former  had  a size  range  from  O.Olym  to  lO.Oym  and  the  latter 
a size  range  from  0.5pm  to  25.0pm,  Yamamoto  used  an  averaged  value  of 
100  cm"^  per  level  in  number  concentration  or  0.063  g m“3  in  liquid  water 
content,  whereas  Platt  had  a value  of  only  40  cm-^  or  a liquid  content  of 
0,035  g m“3. 


SUMMARY  AND  CONCLUSIONS 

In  this  study,  a detailed  investigation  was  performed  in  the  10.5pm  to 
12.5pm  window  region  of  the  internal  structure  of  the  vertical  upward 
radiation  emerging  from  the  different  levels  inside  a realistic  inhomoge- 
neous cloud  in  the  hope  to  derive  additional  physical  insight  for  an  under- 
standing of  cloud  radiation.  The  contributions  to  the  internal  radiation 
field  consists  of  two  components:  one  includes  the  thermal  emission  of 
cloud  particles  and  water  vapor  and  the  resulting  Mie  scattering  by  these 
particles,  and  the  other  includes  the  thermal  emission  from  the  ground  or 
sea  surface  and  the  resulting  scattering.  As  the  former  increases  with 
cloud  depth,  the  latter  decreases.  When  both  are  present,  the  cloud  is 
not  "black"  and  its  radiation  field  is  said  to  be  distorted  by  the  pene- 
tration of  surface  radiation  and  scattering.  Such  penetration  is  greater 
at  10.5pm  than  at  longer  wavelengths.  When  the  penetration  ceases,  the 
cloud  emission  and  scattering  completely  dominate  the  radiation  field  and 
the  cloud  then  behaves  like  a blackbody  radiator.  The  fractional  amounts 
of  thermal  emission  and  scattering  contributing  to  the  blackbody  radiation 
at  the  cloud  top  are  determined  by  the  average  single  scattering  albedo 
of  the  cloud. 

The  average  meteorological  and  microphysical  properties  were  utilized  to 
examine  clouds  of  various  homogeneous  states  at  two  wavelengths  for  their 
internal  radiation  fields.  The  microphysical  property  of  clouds  was  shown 
to  be  of  prime  importance  in  cloud  modeling,  whereas  the  temperature  dis- 
tribution within  the  cloud  had  little  effect  on  either  the  upwelling  radi- 
ation from  the  cloud  top  or  the  internal  radiation  distribution. 

From  such  detailed  investigation,  a technique  was  evolved  which  takes 
advantage  of  the  distorted  radiation  field  in  a cloud.  A 10.5pm  channel 
is  proposed  in  addition  to  the  existing  window  channel  in  the  12pm  window 
region  in  satellite  sounders.  The  ratio  of  the  observed  upward  radiances 
at  these  two  wavelengths,  when  compared  with  that  of  the  blackbody  radiances 
derived  for  the  corresponding  wavelengths  from  a knowledge  of  the  surface 
temperature,  which  could  be  an  estimate,  offers  a good  indication  of  the 
degree  of  the  cloud  blackness  or  of  cloud  emissivity.  Consequently,  it 
appears  possible  to  determine  the  reliability  of  the  temperature  field 
derived  from  satellite  data.  It  was  shown  that  the  technique  for  distin- 
guishing the  cloud  blackness  does  not  require  an  accurate  surface  tempera- 
ture. Two  examples  were  presented  in  which  one  cloud  was  not  so  "black" 
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and  the  other  was.  The  clouds  which  have  satisfied  the  ratio  test  at 
different  levels  normally  have  emissivity  values  about  0.95  or  more.  On 
the  other  hand,  when  the  observed  upward  radiances,  upon  conversion  to 
temperature,  gtve  nearly  the  same  value,  cloud  emissivity  may  be  taken 
to  be  unity  and  that  temperature  then  the  cloud-top  temperature. 

Finally,  this  proposed  technique  is  simple  and  inexpensive  to  execute  on 
board  a satellite.  Moreover,  the  ratio  test  minimizes  the  systematic 
errors,  which  are  found  in  all  measuring  devices  and  which  are  difficult 
to  account  for. 
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